Quantum phases in a resonantly-interacting Bose-Fermi mixture 
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We consider a resonantly-interacting Bose-Fermi mixture of 40 K and 87 Rb atoms in an optical 
lattice. We show that by using a red-detuned optical lattice the mixture can be accurately described 
by a generalized Hubbard model for 40 K and 87 Rb atoms, and 40 K- 87 Rb molecules. The microscopic 
parameters of this model are fully determined by the details of the optical lattice and the interspecies 
Feshbach resonance in the absence of the lattice. We predict a quantum phase transition to occur 
in this system already at low atomic filling fraction, and present the phase diagram as a function of 
the temperature and the applied magnetic field. 

PACS numbers: 03.75.Fi, 67.40.-w, 32.80.Pj, 39.25+k 



Introduction. — In the last few years experiments have 
shown that it is possible to realize a quantum degenerate 
gas of fermionic atoms p| . Combining such a degenerate 
Fermi gas with a Bose-Einstein condensate of atoms, it is 
also possible to obtain a quantum degenerate Bose-Fermi 
mixture Q, thus creating a dilute analog of the liquid 
3 He- He mixture. The recent observation of interspecies 
Feshbach resonances in a Bose-Fermi mixture 0, opens 
up even richer physics, as this couples the fermionic and 
bosonic atoms to a third species, namely the fermionic 
heteronuclear molecules. An interesting aspect of the 
resonantly-interacting mixture that we address in much 
more detail in the following, is the possibility to re- 
versibly destroy a Bose-Einstein condensate of atoms, by 
associating the bosonic atoms into fermionic molecules 
and thus creating a degenerate Fermi gas of dipolar par- 
ticles. Moreover, the Bose-Einstein condensed phase is 
very interesting by itself because it contains a macro- 
scopic quantum coherence between the fermionic atoms 
and molecules. However, in experiments with magnetic 
or optical traps the molecules quickly decay due to atom- 
molecule and molecule- molecule collision |5(. By loading 
the degenerate Bose-Fermi mixture into an optical lattice 
with a total filling factor less then unity, these collisions 
can be prevented and the lifetime of the molecules is ex- 
pected to be dramatically enhanced. 

Bose-Fermi mixtures in an optical lattice, but in the 
absence of an interspecies Feshbach resonance, have been 
the subject of active theoretical investigation lately. In 
particular, lattice symmetry breaking 6], the existence 
of quantum phases that involve the pairing of fermions 
with bosons 0], and also unconventional fermion pair- 
ing [j| have been predicted. It is the main objective 
of this Letter to generalize these studies to the case of 
a resonantly-interacting Bose-Fermi mixture. Although 
our theoretical methods are very general, we consider as 
a concrete example a mixture of fermionic 40 K atoms and 
the bosonic 87 Rb atoms, as this system is now becoming 



available experimentally 9]. A gas consisting of these 
two atoms is especially promising as their wavelength is 
relatively close and readily accessible using Ti:Sapphirc 
lasers. Furthermore, the mass of K is much larger than 
the mass of the other experimentally available fermionic 
atom 6 Li, which makes it much easier to trap this species 
in an optical lattice. 

In order to analyze the properties of a resonantly- 
interacting Bose-Fermi mixture most easily, it is conve- 
nient that all the species in the optical lattice, i.e., the 
fermionic atoms, the bosonic atoms, and the fermionic 
molecules, experience the same on-site trapping fre- 
quency. Because the potassium atoms are lighter than 
the rubidium atoms, the potential for the former should 
thus be less deep than that for the latter. Moreover, 
as is shown in the inset of Fig. ^ both the Dl and D2 
lines of potassium are blue compared to the Dl and D2 
lines of rubidium. As a result equal on-site trapping fre- 
quencies can only be achieved in an optical lattice that 
is red detuned with respect to all these four transitions. 
Making use of the fact that the hyperfine structure of 
the atoms is no longer resolved for the detunings used in 
optical lattice experiments we have calculated the 
ratio of the two on-site trapping frequencies as a func- 
tion of the wavelength of the lattice laser. As can be 
seen from Fig. ^ using a wavelength of 806 nm, ensures 
that the trapping frequencies for both atomic species are 
the same. In principle the polarizability of the molecule 
is not known. However, recent experiments have shown 
that for homonuclcar molecules, the resulting trap fre- 
quency for the molecules is almost the same as that of 
the atoms [llj . In the following, we make the reasonable 
assumption that this also holds for the heteronuclear 40 K- 
87 Rb molecules of interest to us. 

For the particular detuning given above, the on-site 
trapping frequency is related to the Rabi frequency O of 
the lattice laser by u> = ujf = = 2.1 • 10~ 5 Q. Having 
to use a red-detuned optical lattice has the disadvantage 



2 



that the atoms are trapped in the light and not in the 
dark, which in principle results in a larger decay due to 
spontaneous emission. We have therefore also calculated 
this emission rate and find that T /VI 2 — 1.9 • 1CU 21 s for 
the optimal wavelength. Clearly the lifetime is always 
longer than one second for realistic Rabi frequencies less 
than 4 GHz, and spontaneous emission can thus be safely 
neglected. 
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FIG. 1: Ratio of the on-site trap frequencies of the boson 
87 Rb and the fermion 40 K, as a function of the lattice laser 
wavelength. The inset shows the fine-structure levels of the 
two atomic species. 

Generalized Hubbard model. — We now consider the 
many-body aspects of the Bose- Fermi mixture near a Fes- 
hbach resonance. We consider a mixture of 40 K and 87 Rb 
that is loaded into a three-dimensional and cubic optical 
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lattice that is sufficiently deep, so that we are allowed 
to use a tight-binding approximation for the band struc- 
ture of the single-particle states. Experimentally this re- 
quires the Rabi frequency of the lattice laser to be larger 
than about 1 GHz. For the reasons mentioned above, 
we consider only low filling fractions, which limits the 
maximum number of atoms on a single site. As a result, 
we can neglect atom-molecule and molecule-molecule in- 
teractions. Furthermore, for low filling fractions there is 
also no problem of phase separation of the atomic Bose- 
Fermi mixture |l2| and it is just ified to neglect possible 
Mott physics in the mixture |13l 1 1 - JJ ] . 

Under these conditions we have recently derived the 
theory for resonantly-interacting ultracold atomic gases 
in an optical lattice 0]. Within this theory the two- 
body Feshbach problem at a single site is solved exactly, 
which physically leads to a dressing of the molecules and 
to various avoided crossings in the on-site energy levels of 
two atoms e a (B) as the magnetic field B is swept through 
the Feshbach resonance. After having solved the on-site 
problem the various hopping parameters can be calcu- 
lated in the tight-binding approximation. In this man- 
ner the microscopic parameters of the generalized Hub- 
bard model that describes the Bose-Fermi mixture near 
the Feshbach resonance are completely determined by the 
details of the optical lattice potential and the experimen- 
tally known parameters of the Feshbach resonance in the 
absence of the optical lattice. Ultimately, the mixture is 
described by the following effective hamiltonian, 
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Here tp, is and t m are the tunneling or hopping ampli- 
tudes for the fermionic atoms, the bosonic atoms, and 
the fermionic molecules, respectively. The symbol 
denotes a sum over nearest neighbors. The operators c], 
c i and a\ , a, t correspond to the creation and annihilation 
operators of a single fermionic and bosonic atom at site 
i, respectively. The operators b\ al and b ia correspond 
to the creation and annihilation operators of the vari- 
ous dressed molecules at site i. Also e a = 3hu>/2 is the 
on-site energy of a single atom. The on-site interaction 
between two bosons is given by U BB , and U BF denotes 
the on-site background interaction between the bosons 
and the fermions. In the tight-binding limit the hopping 
amplitudes can be conveniently expressed in terms of the 



lattice parameters as 

Here v distinguishes the different species in the mixture, 
i.e., v = F for the fermionic atoms, v = B for the bosonic 
atoms, and v — m for the molecules. Moreover, A is the 
wavelength of the light used to create the optical lattice, 
and the harmonic oscillator lengths obey l v = -^/h/m^uj, 
where mp, me, and m m = mp + tob are the masses 
of the fermions, bosons, and molecules, respectively. As 
expected, we have in general that £ m <C £f,b -C tuo. 
Note that we have introduced a chemical potential for 
each atomic species, since it is experimentally possible to 
control both the filling fraction of the fermions as well as 



3 



o - 



-4 - 




-0.05 0.05 
B - B (Gauss) 

FIG. 2: Details of the physical content of our theory. We 
show the avoided crossing between the bare molecular level 
and the lowest two-atom trap state. The inset shows the 
probability Z a as a function of the magnetic field. This figure 
was calculated for the 40 K - 87 Rb mixture with A = 806 nm 
and Q = 1 GHz. 



the bosons in the mixture. 

Close to resonance we can neglect the on-site interac- 
tions. The strength of the atom-molecule coupling in 
the lattice is given by g' = g/(Tr(l B + ^)) 3 ^ 4 , where 
g = ^/27rabgA/iAi?/r7i r is the bare atom-molecule cou- 
pling without the lattice and m r = uiFmB / '(tuf + tub) 
is the reduced mass. For the 40 K- 87 Rb mixture with 
potassium in the hyperfme state 1 9/2 , —9/2) and rubid- 
ium in the hyperfme state |1, 1) there occurs a Feshbach 
resonance at Bq = 510 Gauss for which the parameters 
that determine g are given by the background scatter- 
ing length dbg = 150 ao and the width of the resonance 
AB = 1 Gauss 0, 0] . The difference in magnetic mo- 
ments Afi is equal to 29/22 Bohr magneton in that case. 
In Fig. |^1 we show a close-up of the avoided crossing and 
the wavefunction renormalisation factors Z a that give the 
probability for the dressed molecules to be in the bare 
molecular state of this Feshbach resonance. The proba- 
bility Z a is determined by the selfenergy of the molecules 
HE m (E) = (g 2 /i:lfhuj) ■ T{-E/2fiuo + 3/4)/T{-E/2huj + 
1/4), with r(z) the gamma function, through the rela- 
tion Z a = 1/(1 -dhE m (E)/dE) [13. Note that in Fig.H 
the sum of the probabilities + Z^ does not add up 
to one. This means that for this relatively broad inter- 
species Feshbach resonance a single-band approximation 
is not valid to determine the dressed molecular wavefunc- 
tions. However, a single-band approximation in terms of 
dressed molecules is always possible for the low filling 
fractions of interest to us, because the higher-lying on- 
site dressed molecular states will not be populated as we 
show now. 

Phase diagram. — With the above formalism we next 
determine the phase diagram for the K and Rb mix- 



ture. For simplicity we consider only equal densities for 
both atomic species, although the generalization is imme- 
diate. We mentioned already that we consider a deep op- 
tical lattice for which the hopping strengths of the atoms 
are small with respect to the level splitting in the on-site 
microtrap realized by the optical lattice. The hopping 
strength of the molecules can be completely neglected in 
this limit and consequently the band-structure for the 
molecules is essentially flat. At zero temperature the en- 
ergy of the lowest molecular state can, depending on 
the magnetic field, be either smaller or larger than the 
sum of the lowest atomic energy levels. As a result, the 
ground state is either a Fermi sea of 40 K - 87 Rb molecules, 
or a Fermi sea of 40 K atoms and a Bose- Einstein conden- 
sate of 87 Rb atoms. Because of the different symmetries 
of these ground states there exists a quantum phase tran- 
sition between these two states that is in the same uni- 
versality class as the XY model. In the following, we 
calculate the phase diagram as a function of total fill- 
ing fraction and temperature by performing a mean-field 
analysis of the normal state of the hamiltonian in Eq. 
©■ 

For the equation of state for the total filling fraction 
we find always that 

n = n F + n B + 2^n m!(T . (3) 

a 

In the normal state the molecular filling fractions obey 



= J_V l - 
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where Titu^a — e^ a — (Mf + ^b) is the molecular disper- 
sion. The atomic filling fractions are given by 
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where ?jWk,F = £k,F — and Tiwk.B = £k.B — fJ-B 

are the dispersion relations for the fermionic 
and the bosonic atoms, respectively. Here 
£k,F,B = — 2ip,B J2j=i cos (%V2) + £ a for the atoms 
and ek.o- = for the molecules. We have seen that at 
zero temperature there is a quantum phase transition 
between a phase where the 87 Rb atoms are Bose- 
Einstein condensed and a phase with only a Fermi sea of 
molecules. The quantum critical point is determined by 
the ideal gas condition for Bose-Einstein condensation, 
i.e., e| = 2e a — 6(£b + £f)- hi the approximation that 
we can neglect the hopping of the molecules, their band 
structure is flat and the quantum critical point is inde- 
pendent of the filling fraction of the molecules. Including 
the molecular band structure would lead to a critical 
magnetic field that slowly shifts to lower magnetic fields 
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FIG. 3: Phase diagram of the 40 K and 87 Rb mixture as a 
function of the total filling fraction, magnetic field and tem- 



perature. This figure was calculated for the K 
ture with A = 806 nm and Q — 1 GHz. 



Rb mix- 



magnetic field and temperature. Our analysis of the 
quantum phases of a resonantly-interacting Bose-Fermi 
mixture has been based on mean-field theory. In partic- 
ular, this implies that we have not considered the attrac- 
tive finite-range interaction between the fermionic atoms 
that can be mediated by density fluctuations in the Bose- 
Einstein condensate In principle this mechanism 

can lead to a BCS pairing between the fermionic atoms. 
However, for the spin-polarized mixture discussed here 
this pairing must take place in a p-wave channel, which 
is expected to have a very small critical temperature at 
small filling fractions. We, therefore, do not consider this 
interesting possibility in detail here and leave that to fu- 
ture investigations. 

This work is supported by the Stichting voor Fun- 
damenteel Onderzoek der Materie (FOM) and by the 
Nederlandse Organisatie voor Wetenschaplijk Onderzoek 
(NWO). D. van Oosten acknowledges funding through a 
Marie-Curie Excellence grant. 



as the density increases. At nonzero temperatures we 
can also determine the critical temperature as a function 
of the detuning and the total filling fraction from the 
equation of state. The phase diagram in Fig. [21 shows 
how at constant total atomic filling fraction the critical 
temperature depends on the magnetic field. For large 
enough magnetic fields there are no molecules and the 
critical temperature is determined for low densities 
by the critical temperature of an ideal gas, which is 
proportional to n 2 ^ 3 . Note that the critical temperature 
always obeys T c <C Tiw/fcs, which a posteriori shows that 
a single-band approximation for the dressed molecules is 
indeed appropriate. 

Conclusions and Discussion. — In summary we have 
shown that by using a red-detuned optical lattice with 
a wavelength of 806 nm the Bose-Fermi mixture of 40 K 
and 87 Rb atoms can be accurately described by a gen- 
eralized Hubbard model. Moreover, we have shown that 
the model contains a quantum phase transition associ- 
ated with the Bose-Einstein condensation of rubidium. 
Interestingly, the presence of a Bose-Einstein conden- 
sate induces also a macroscopic coherence between the 
fermionic atoms and molecules, because of the specific 
form of the atom-molecule coupling near a Feshbach res- 
onance. What is especially interesting is that such a co- 
herence cannot be obtained by solely making use of lasers 
in this case because it involves a quantum coherence be- 
tween two different species. It would, therefore, be very 
exciting to observe Rabi oscillations between fermionic 
atoms and molecules by an appropriate manipulation of 
the atomic Bose-Einstein condensate density. 

Using the known atomic physics of the Feshbach res- 
onance to determine the parameters in the generalized 
Hubbard hamiltonian we have calculated the critical sur- 
face for low filling fractions as a function of the applied 
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